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Abstract. The contribution to the ultraviolet background (UVB) from thermal 
emission due to gas shock heated during cosmic structure formation is assessed 
with an updated version of Press-Schechter (Sheth & Tormen 1999) formalism. 
The calculation is consistent with empirical estimates based on the observed prop- 
erties of galaxies and the observed cosmic star formation history. The bulk of the 
radiation turns out to be produced by objects in the mass range lO'^^"'^^ A^q , i.e. 
large galaxies and small groups. When compared to more conventional compo- 
nents (QSOs, stellar) it is found that near 1 Ry, thermal emission is comparable 
to stellar contributions and amounts to about 10 %, 20 % and 35 % of the total 
flux at redshifts of 3, 4.5 and higher, respectively; more importantly, near the ion- 
ization threshold for He ii , the thermal contribution is comparable to the QSO 
intensity already at redshift ~ 3 and dominates at redshifts above 4. In addi- 
tion, thermal photons alone are enough to produce and sustain He ii reionization 
already at 2 ~ 6. Finally, the observed Gunn-Peterson effect at high redshifts 
(3 < z < 6) constrains the escape fraction of ionizing photons from galaxies to 
less than a few percent. 
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1. Introduction 

Neutral hydrogen in the intergalactic 
medium (IGM) produces a forest of reso- 
nant Lya absorption lines in the spectra 
of high-redshift quasars. The connection of 
these features to the structure formation 
process has been now firmly established us- 
ing N-body/hydrodynamic n umerical simu- 
lations ifCen etrar.lii99 4: Zhan g et alJll99nl: 
iMiralda-Escude et al.l>199fi[l . There is con- 
sensus that the observed IGM tempera- 
ture results from a balance between photo- 
ionization heating and adiabatic cooling 



due to the Hubble expansion. Such photo- 
heating is provided by the extragalactic ul- 
traviolet background (UVB) whose nature, 
origin and evolution have, therefore, been 
subject to considerable investigation. 

At low redshifts, the ionization balance 
is consistent with a pure power-law UVB. 
Traditionally QSOs have been considered 
as the main sources of ionizing photon. 
However, a nu mber of recent resu lts are at 
odds with this. iKim et all l|2n0lD find the 
break in the redshift evolution of absorbers 
at lower redshifts than predicted by numer- 
ical simulations using a standard QSO ion- 
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izing background l|Haardt fc MadaiJUoOfih . 
hinting at an in complete descr i ption o f the 
UVB. This led iBianchi et af] EcM BOl 

hereafter) to recompute the UVB as a 
superposition of contributions from QSOs 
and galaxies. 

Similarly, high resolution simulations 
using a cold dark matter model and a 
standard QSOs ionizing background, pro- 
duce Ly-a forest lines with a minimum 
width signific a ntly b elow that observe d 
^Theuns et a,lJ ^M, iBrvau et a,lJ ^^M- 
This has prompted several suggestions for 
additional heat sources, including photo- 
electric dust heating (Nath et al. 1999), ra- 
diativ e transfer effects (fAbel fc Haehneltl 
Il999fl . or Compton heating by X-ray 
backg round photons l|Madau fc EfstathioiJ 

The situation around redshift z ~ 
3 is more complicated. There is a gen- 
eral trend for the optical depth at the 
He"*" edge to inc r ease with redsh ift (e.g. 
iHeap et al.ll200nD . ISongailal l|l99a^ reports 
an abrupt change of the C iv /Si iv ratio 
at z PS 3, often interpreted as a hardening 
of the ionization spectrum due to a sud- 
den He II reionization. This is corroborated 
by the detection of patchy He ii Ly—a 
absorption at similar re dshifts. However , 
more recent VLT/ UVES (iKim et a l. 2002) 
and Keck/HIRES l|Boksenberg et al 2003) 
studies find no such discontinuity around 
z — 3, and even suggest that C iv /Si iv is 
not a good indicator of the He ii ionization 
state. 

Here we report on recent findings 
about another source of UVB ionizing pho- 
tons, namely thermal emission from shock- 
heated gas in collapsed cosmic structures. 
We show that the ionizing photons emit- 
ted by this process make a non negligi- 
ble fraction of the metagalactic flux, that 
the resulting spectrum is hard and copious 
He II ionizing photons are produced. In fact 
this process may well dominate the produc- 
tion of such hard photons at z > 4. 



2. Model 

We consider the mean ultraviolet ioniz- 
ing flux produced by QSO, stellar and 
thermal components. After estimating the 
mean volume emissivity of each component 
as a function of redshift, as described be- 
low, we solve for the cosmological radiative 
transfer equation (Peebles 1993). The ef- 
fective optical depth, t'^^ , due to the pres- 
ence of intervening neutral gas, is com- 
puted by assuming a distribution of ab- 
sorbers as a f unction of H i colum n density 
and redshift ijParesce et al.lll985() . The de- 
pendence on Nm is determined from counts 
of Lya absorption lines in QSOs spectra; 
a power-law evolution with redshift such 
that r°*^ (1 -I- zy, is further assumed 
(jin^lilil- We use 7 ~ 3.4 in the red- 
shift range 1.5 < z < 4 l)Kim et al.ll200H) . 
and, based on the strong evolution for r®*^ 
as impl ied by recent spectra of high redshift 
QSOs iBecker et al.ll200lt iFan et al.ll20(il 
we adopt 7 = 5.5 for z > 4. 

The contributions of QSOs and stars 
to the ionizing UV background adopted 
here are similar to those in BOl. Basically 
the QSO emissivity assumes a luminosity 
function that follows the double power-law 
model of Bovlc et al. (1988). For z < 3, we 
adopt the parameters given in iBovle et aTl 
1)2000(1 . At z > 3, we inc lude the exponen- 
tial decline suggested bv lFan et all JloOl). 
For the QSO spectrum in the ionizing UV 
range we use a sir nple p ower law, j{u) oc 
ly-^-^ llZheng et aH[l997i) . 

For the stellar component, we as- 
sume a star formation rate that is con- 
stant from high redshifts to z w 1, 
then rapidly decreases to local values, 
as indicated by galaxy su rveys in the 
rest-frame non-ion izing UV l|Madaii et al.l 
ll99aHSteide l et al."l99;j^. Synthetic galac- 
tic spectra (Bruzual A. & Chariot 1993; 
version 2001) are then used to calculate the 
emissivity of ionizing UV photons as a func- 
tion of z. The absorption of radiation by 
the galaxy interstellar medium is modeled 
by a redshift-independent value for /esc, the 
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fraction of Ly-continuum photons that can 
escape into the IGM. 

Finally, for the emissivity from shocked 
intergalactic gas we assume such gas to be 
thermalized in virialized dark matter halos, 
which we de scribe with a Press-Sch echter 
prescription l|Sheth fc Tormenl ll999D . The 
dark matter density in virialized objects is 
computed according to the spherical col- 
lapse model; the ratio of baryons to dark 
matter is taken as the universal value, ~ 
0.15. The spectral thermal emissivity of the 
shocked gas in each ha lo is thus computed 
throug h the code by iRavmond fc SmitiD 
l)1977l version 1992; collisional regime is a 
good approximation for the plasma condi- 
tions); we then estimate the thermal vol- 
ume emissivity by adding up the contribu- 
tions from all collapsed structures. We ac- 
count for the finite cooling time of a halo, 
by introducing a term ensuring that the ra- 
diated energy does not exceed the thermal 
energy of the system. Finally, we take into 
account feedback effects by computing the 
increase in the cooling time caused by the 
injection of additional, feedback energy. 

More de ta ils are presented in 
iMiniati et all lj2004|) . We note that 
cooling process associated with thermal 
emission investigated here also leads to 
galax y /star f orination__j^^i^^^^ei 
11978(1 ■ In fact, iMiniati et alJ l|2004ri also 
show that the estimated flux presented 
here, based on a Press-Schechter model, is 
in remarkable agreement with an empirical 
estimate based on the observed cosmic 
star formation history and the distribution 
of stellar mass as a function of halo virial 
temperat ure, as reconstructed fro m SDSS 
data from lKauffmann et al.l l|2003j) . 

3. Results 

3.1. Photo-ionization Rates 

Fig. n] shows the evolution as a function 
of redshift, of the photo-ionization rates in 
units of 10~^^ s~^ defined as 



for the various UV radiation compo- 
nents discussed above, together with mea- 
sured values inferred from the observed 
Gunn-Pctcrson effect in hig h-z QSO spec- 
tra (.Mcpona,1d Mira,lda,-Escudel 120011 
IFan et al.ll20m 

The left panel of Fig.^is relative to ion- 
ization of H I and contains a number of im- 
portant features. First, although the emis- 
sion from QSOs is able to produce the ion- 
izing flux observed at 2 — 3, it falls 
short at higher redshifts, a well known fact. 
In our formulation, it results from the as- 
sumed rapid decline of the QSO number 
dens ity for z > 3, as d erived from the SDSS 
(33 IFan et al.ll200l . 

Second, a comparison of the measured 
values of and the stellar ionization 
rates assuming fesc — 1% (dash curve) im- 
plies, according to our model, that fesc is 
smaller than a few %. Estimates of the UV 
background from the proximity effect are 
known to be larger than those obtained 
via theoretical mod els of the IGM opac- 
ity (i.e. the work of iFan et alJ we are us- 
ing here), probably because of a bias of the 
QSO distri bution towards the dense r envi- 
ronments ((Schirber fc Bullockll2002(l or be- 
caus e of systema tic errors due to line blend- 
ing l|Scott et al .11200(1 . 

Finally, the right panel of Fig. ^ shows 
the He ii ionization rates. According to 
the plot, above 4 Ry stellar emission is 
thoroughly negligible (independent of /esc) 
whereas thermal emission is comparable to 
QSOs at z ~ 3 but completely dominates 
the radiation flux at higher redshifts. This 
result is very important in terms of the 
IGM evolution and has not previously been 
noticed. It depends only weakly on feed- 
back, but it does assume an escape fraction 
of the thermal photons from collapsed ha- 
los of order 1. As discussed in Miniati ct al] 
l)2004(l these conditions should be ensured 
for coUisionally ionized gas within halos of 
virial temperature above 10^ K, that is for 
the halos that generate most of the thermal 
emission. 



26 



F. Miniati et al.: UV Flux from Structure formation 




10- 



1 0- 



1 





1 1 1 1 1 1 1 1 

Hell _ 




\^^^ — 
\ ^\ ; 
\ : 


xlO 


\ ^^^^ 






1 


\ \ 

1 1 1 1 1 1 1 1 \i 1 


3 


4 5 6 
z 



Fig. 1. Left: Photo- ionization rates defined in eq-^as a function of redshift for ionization 
of H I due to emission from: QSO (dash), stellar for f^sc — 1% (long dash) and shocked 
IGM with feedback effects (solid) and without (dot). The thick solid line is the total 
considering the feedback case (or, alternatively, the no-feedback case and fesc = 2%). 
The data points are from[McDonald_&_Mir.alda-Escudc (2001) and Fan et.al. (20J22,) after 
correction for our cosmological model. Right: same as for left panel but now for coefficients 
relative to He ii (right). The stellar component (long dash) has now been multiplied by 
a factor 10 for visualization purposes. 



3.2. Softness parameter 

Both for the individual spectral compo- 
nents and for the total spectrum, we have 
computed the redshift evolution of the 
spectral softness parameter, Sl = -j^i 
where Jhi and Jneii are the UVB in- 
tensities at the H i and He ii Lyman lim- 
its respectively. These are shown in Fig. 
121 Thermal emission are characterized by 
moderate values of Sl, raising with red- 
shift from ~ 10 to 100 and 80 to 300, 
respectively; the stellar component shows 
less evolution but with maximal Sl val- 
ues above 1000 at z ^ 3. The composite 
spectra (feedback case) Sl evolution re- 
sembles very closely that of QSOs, even 
at high redshift: this is somewhat fortu- 
itous as at 2; ^ 4 the UVB is dominated 
by the sum of stellar and thermal contri- 
bution. When compared with t he available 
data, e.g. the recent study of iHeap et al.l 
II2000D . our predicted values for S'l at red- 



shift around z = 3.2 are close to the mea- 
surements of those authors. However, the 
reported change in the parameter S would 
not necessarily imply a "jump" of the same 
quantity at z ~ 3, as it would be due to 
the evolution of our composite spectrum; 
in particular to the decrease in the ratio of 
thermal to the QSO flux. 

3.3. He II Reionization 

The previous results hint at the intrigu- 
ing possibility that He ii reionization could 
have been powered by UV light from cos- 
mic structure formation. For this to be the 
case, the production rate of ionizing pho- 
tons has to satisfy the condition 

F X min(Trec,T//„6f,;e) > 1. (2) 

Here Tree — 0-Q/ngas<x{T)C is the recom- 
bination time, a is the radiative recombi- 
nation coefficient, and the clumping factor 
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Fig. 2. Redshift evolution of the softness 
ratio Sl for different emission cases: ther- 
mal, no feedback (dot); thermal, with feed- 
back (thin solid); stellar (long dash); QSOs 
(short dash); total, no feedback (higher 
thick solid); total, with feedback (lower 
thick solid) 

C = {n^) / {np)"^ > 1 is meant to allow for 
the effects of density inhomogeneities in- 
side the ionized region. Using a helium to 
hydrogen number ratio y = 0.08 and as- 
suming a temperature of the reionized gas 
T « 4 X 10'' K, we find 

Tree = 5.3 X lO^^ C'^ 

which is shorter than a Hubble time for 
z ^ 4.5. Thus, from Fig. ^ we find that at 
z « 6 thermal emission dominates the pho- 
toionization rates and alone provides 5.5 x 
10~^^ (10~^^) He II photoionizations/s in 
the no-feedback (feedback) case. According 
to eq. (PJ, the He ii recombination rate 
at the same redshift is 6.4 x 10^ ^^C s~^. 
Hence, it appears that structure forma- 
tion can produce He ii reionization around 
z = 6, without the contribution from any 
other process and essentially independently 
of the feedback prescription adopted. 



Whether or not this possibility is fully 
compatible with all observational results is 
not clear at the moment and will be the 
subject of future investigations. 

4. Summary 

We have shown that UVB ionizing pho- 
tons can be copiously produced by thermal 
emission from shock-heated gas in collaps- 
ing cosmic structures. 

Thermal radiation is characterized by 
a hard spectrum extending up to pho- 
ton energies of order hpV ^ ksT. This 
is well above the H i and He ii ionization 
thresholds for virial temperatures above 
10^ K. The bulk of the emission is pro- 
duced by halos with temperatures between 
10^ K and a few xlO'^ K, corresponding 
to masses 10^^""'^^ Mq. We assume that 
most of the thermal radiation is able to 
freely escape into intergalactic space, which 
is justified for a gas that is coUisionally ion- 
ized and at th e temperature of these halos 
l|Miniati et al.l [20041 . 

We use simplified radiative transfer to 
compute the transmitted flux due to QSO, 
stellar and thermal emissions. Importantly, 
the resulting associated photoionization 
rates, when compared to measurements 
of the Lyman series Gunn-Peterson ef- 
fect in the spectra of high redshif t QSOs 
l|Fan et al.ll2003lBecker et al.ll200ll) . imply 
an escape fraction of UV ionizing photons 
from galaxies, /esc, below a few %. This re- 
sult is in agreement with very recent and 
indepe ndent determinations o f /p.^r, c arried 
out by iFernandez-Soto et alJ <)2003j) , who 
set a 3(7 (statistical) upper hmit fesc ~ 4% 
for a sample of spectroscopically identified 
galaxies of redshift 1.9 < z < 3.5 in the 
Hubble Deep Field. 

With fesc — 1%, it turns out that near 
the H I ionization threshold, thermal emis- 
sion is comparable to the stellar compo- 
nent and amounts to about 5-10 %, 15- 
30 % and 20-50 % of the total at redshifts 
of 3, 4.5 and higher respectively. Near the 
ionization threshold for He ii , the ther- 
mal contribution is much stronger. It is 
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comparable to the QSO input already at 
2 ~ 3, and it dominates for z > 4. Thus, 
this contribution, with a typical softness 
parameter Sl = 10 — 100, is expected 
to play a major role in He ii reionization. 
In principle structure formation alone pro- 
vides enough photons to produce and sus- 
tain He II reionization at 2 ~ 6. 
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